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Abstract — In this letter, the problem of implementing compute 
and forward (CF) is addressed. We present a practical signal 
model to implement CF which is built on the basis of Gaus- 
sian integer lattice partitions. We provide practical decoding 
functions at both relay and destination nodes thereby providing 
a framework for complete analysis of CF. Our main result is 
the analytical derivation and simulations based validation of 
union bound of probability of error for end to end performance 
of CF. We show that the performance is not limited by the 
linear combination decoding at the relay but by the full rank 
requirement of the coefficient matrix at the destination. 

Index Terms — Compute and Forward, Gaussian integers, finite 
fields. 



I. Introduction 

In wireless networks with multiple users, relaying is an im- 
portant technique adopted to maximize the network through- 
put. In HI, Nazer and Gastpar proposed a novel strategy 
of generalized relaying called Compute and Forward (CF) 
which enables the relays in any Gaussian wireless network 
to decode linear equations of the transmitted symbols with 
finite field coefficients, using the noisy linear combinations 
provided by the channel. The linear equations in finite field 
are transmitted to the destination and upon receiving sufficient 
linear equations, the destination can decode desired symbols. 
Further, information theoretical tools are used in JT] to ob- 
tain the achievable rate regions. An algebraic approach to 
implement CF has been introduced in [2| where the authors 
propose to implement CF making a connection between CF 
and isomorphism in module theory. 

The main contribution of this correspondence is to demon- 
strate the implementation of CF using practical signal con- 
stellations and study its end to end performance from source 
to destination. We use signal constellations based on one 
dimensional Gaussian integer lattices to implement CF. We 
utilize the natural isomorphism existing between these signal 
constellations and finite fields (|3l. 0) and apply it to design 
practical encoding and decoding functions at each node of 
the system from source to destination. In order to understand 
the factors affecting the CF behavior, we consider integral 
channels. Therefore, we bypass the errors introduced due to 
non-integral nature of the channel thereby avoiding the "self- 
noise" [ 1 1. We show that at high SNR, the overall performance 
of CF is determined primarily by the choice of the finite field 
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used and is not limited by the detection of linear combinations 
at the relay. We also provide a tight union bound estimate of 
probability of error at the destination of CF. 

II. Preliminaries : Gaussian Integers 

In this section, we will present some useful algebraic 
preliminaries relevant to this letter. Details can be found in 
0, 0. 

Let Q be the Gaussian Integers Z[i] and let Q n denote the 
residue class Q modulo ir where 7r e Q. Any element of Q 
can be mapped to the residue class using the function pL : 
Q — » Q n . which is defined as 



Ms) = 9 



g.ir 



(1) 



where ir* is the conjugate of ir, and [.] is the rounding 
operation which is defined on complex numbers as [a + bi] = 
[a] + [b]i. The analogy of Q and Q w in integer domain is Z 
and Z p for some modulo residue class Z mod p. 

The Gaussian primes are the primes in Gaussian integers 
which are given by (i) ±1 and ±i, (ii) the rational primes p 
with p = 3 mod 4 and (iii) the factors a + ib of rational primes 
p with p = 1 mod 4. The Gaussian primes of type (iii) exist 
for every p = 1 mod 4 because the rational primes of type 
p = 1 mod 4 can be written as sum of squares a 2 + b 2 by the 
well known Fermat's Theorem [5, Pg. 291]. Therefore, 

p = a 2 + b 2 = (a + bi)(a - bi) 

In this letter, we focus on Gaussian primes of type (iii), 
although extension of this work to other types is straight 
forward. 



III. System Model 

Consider the CF system model with L sources, a relay and a 
destination as shown in figure 1. Let wi € F p be the message 
to be transmitted by the Z-th source (/ = 1,2...L) chosen 
from a finite field F p of order p. The vector of all the source 
messages is given by w = [ wi ... wl ]. Each source 
encodes the message wi into a complex signal constellation 



point using the encoder £ : F p 



to obtain 



xi = £(wi) 



(2) 



The signals are transmitted across the channel to the relay. In 
this model, for the primary understanding, we have assumed 
that the channel gains are Gaussian integers and hence there is 
no "self-noise" due to approximation of channel by an integer 
(T). It is also assumed that channel undergoes slow fading and 
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Fig. 1. End to End CF System Model 



hence remains constant throughout the transmission of each 
signal. The signal obtained at the relay is given by 



y = h\X\ + h 2 x 2 + . . . + h L x L + z 



where hi £ Q is the channel coefficient between transmitter 
I and the relay node, z is i.i.d Gaussian noise given by z ~ 
CJ\f(0,(j 2 ). The signal to noise ratio (SNR) is defined as 

E[\\ *l II 2 ] 



SNR = 



(4) 



The aim of the relay is to compute a linear combination of 
source messages in the original message space v £ ¥ p given 
by 

v = aiwi © a 2 w 2 ■ ■ .aLW L (5) 

where ai £ ¥ p are the linear coefficients chosen on the basis 
of hi and indicates summation over finite field. The estimate 
of v obtained at the relay using the decoder T>r : C — > ¥ p is 
given by 

v = V R {y) (6) 

The estimate of the linear combination v is transmitted to the 
destination. Here we assume this transmission between relay to 
destination is error free and the linear combination is obtained 
at the destination exactly as estimated at the relay. The 
destination obtains L such linear combinations. Therefore, the 
decoder at the destination is given by T>d : {¥ p } L — > {¥ p } L 
such that 

w = V D (v) 

where w is the estimate of the L original source signal 
vector w and v is the vector of estimates of the L linear 
combinations. 

IV. Proposed Encoding and Decoding Functions 

In this section, we propose the encoding function for the 
sources and the decoding functions at the relay and the 
destination in order to implement CF scheme. 

A. Construction of the Signal Constellation 

We define some standard useful functions [3| which we 
utilize in constructing the signal constellations to implement 
CF. A signal constellation feasible to implement CF is desired 
to be isomorphic to a finite field. Therefore, a natural choice is 
the residue class of Gaussian integers Q n because any residue 
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Fig. 2. Constellations from residue class Q n for tt = 2- 
and their respective mapping to finite fields F5 and F13 



- i and tt = 3 + 2i 



(3) class On is isomorphic to a finite field ¥ p if tt is a prime in Q. 



The size of the field is given by p =\ tt | . This isomorphism 



is defined by the bijective function <p : ¥ p 

<f(a) = £ = a ■ 



P 



On defined as 
(7) 



and the inverse ip 1 : Q n — >• F p given by 

a = V" 1 ^) = C-O 71 "*) + TOT!-*) mod p 



(8) 



where u.ir + vir* = 1 and the Euclidean algorithm can be 
applied to calculate u and v. With this isomorphism, Q n and 
F p are mathematically equivalent. In figure 2, some examples 
of residue class Q v along with their finite field mapping are 
shown. We will now propose the encoding and decoding 
functions at the sources, relay and destination. 

B. Encoding at the source 

Let W be the message space which is a finite field compris- 
ing of p elements such that W = ¥ p . The source messages 
are chosen from the message space wi £ W. This message 
space is required to be isomorphic to some complex signal 
constellation S in order to implement CF. The encoding at the 
source is therefore done as follows: 

1. Choose a signal space size as tt = p x l 2 where tt £ Q. 
The signal space is hence given by S = Q n . 

2. For each wi £ W, obtain the isomorphic element in Q n 
using the bijection function in as ip : W —> S such that 

xi = ip(wi) 

The encoded signals are transmitted to the relay where a noisy 
linear combination of the signals is obtained given by (01. In 
the next subsection, we discuss the decoding performed at the 
relay. 



C. Decoding at the relay 

The relay aims to compute the linear combination 

v £ W = F„, 



aiwi © a 2 w 2 . . .a L WL 
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where a; is the finite field mapping of the channel gain hi <^ Q 
given by 

ai = Mhi)) (9) 

Particularly, hi is firstly mapped to the residue class using 
the function /i defined in (1) and then mapped to finite field 
using ip in 0. The decoding process at the relay comprises 
of the following steps: 

1 . From the received signal y, obtain a maximum likelihood 
(ML) estimate of y 



Vail 



are' min II y — t II 
6 tee 11 11 



(10) 



2. Map the ML estimator output with the corresponding 
residue class element in S = Q v using (1) as 



(11) 



The output of this operation yields u 6 Q n which is the 
estimate of linear combination in signal space domain. 

3. Map the estimated signal constellation point u to message 
space given by finite field W = ¥ p using ((HJ to obtain 



(12) 



The output of this operation yields an estimate of the linear 
combination of the original source signals in finite field W. 

An error occurs at the relay if the linear combination is 
incorrectly estimated. More precisely, the probability of error 
at the relay is 

Pr = Pr(u ^ v) (13) 

The relay transmits the estimate of the linear combination to 
the destination where the original source signals are decoded. 

D. Decoding at the destination 

The destination collects L linear combinations from the 
relay which can be written as 



(14) 
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where v l denotes the t— th linear combination (t = 1 . . . L) 
and a\ denotes the Z — th coefficient in t— th linear combination 
between the Ith source and relay given by ©. The decoder at 
the destination inverts the matrix A and obtains an estimate 
of w. Therefore, 

w = A _1 v 

Note that here the inverse of A is taken in ¥ p and A is required 
to be full rank in ¥ p for successful decoding. 

The probability of error at the destination is given by 



V. Probability of Error 

In this section, we derive an analytical expression for 
probability of error at the destination. Since the probability of 
error at the destination is also dependent on the probability of 
error at the relay, therefore, the later is consequently derived. 

Recall from equation ([T5l l that the probability of error at 
the destination is the probability of decoding incorrect original 
source signals such that Pjj = Pr(w ^ w). Therefore, there 
is an error in detection of w, if there is an error at the relay in 
computing any of the L linear combinations of original signals 
or if all the L linear combinations are not independent (and 
consequently, A in (fT~4-b is not full rank). In the next theorem 
, we present a theoretical expression for the union bound on 
the probability of error at the destination. 

Theorem 1. The union bound estimate of probability of error 
at the destination in CF with L sources using finite field 
of size p and Gaussian integer residue class based signal 
constellation is given by 



Pd<Pi + (LPr) 



where 



and 



p,=i-n 



Pr = 1 



erf 



1 



V 2\/2a 



such that a 2 is the variance of additive noise at the relay. 

Proof: An error occurs at the destination if there is 
an error in detection of any linear combination at the relay 
node and/or the linear combinations at the destination are not 
independent (and consequently, A is not full rank). Therefore, 
the union bound estimate of probability of error is given by 

Pd<Pi+^2 P R 

L 

where Pi is the probability of A to have a rank failure (in 
F p ) and Pr is the probability of error at the relay. It has been 
proved in [4] that the probability of an L x L matrix A over 
a finite field of size p, not being full rank is given by 



P 1= Pr(| A|=0) = l-nfl"i 



(16) 



To evaluate the probability of error at the relay, we use the 
classic notion of estimation of error probability. Recall from 
equation ( fT3l that the probability of error at the relay is 
the probability of decoding an incorrect linear combination 
such that Pr = Pr(v ^ w).We rewrite v using dT0b-(TT2b as 
v — </? _1 (^(j/A./i))). Since the maps \i and ip are discrete, the 
equation ( fT3l l can be written as 



P D = Pr(w ^ w) 



(15) 



Pr = Pr(y A1 L ^ (hixi + h 2 x 2 



■ h L x L )) 



Therefore, an error occurs at the destination if the original 
signals are incorrectly estimated. 



Since hi, xi e Q, therefore, the above expression is reduced to 
the probability that the added noise exceeds the voronoi region 
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of Q. The noise is assumed to have a Gaussian distribution with 
mean and variance a 2 . HenceQ 



Pn 



erfc 



( 1 



(17) 



dt. 



where erfc(x) = y= e~ 

Further, the probability of error in decoding L linear com- 
binations at the relay is given by J2l Pr ~ LPr because 
all the transmissions are considered independent. Inserting Pi 
and Pji in union bound estimate, the result is proved. ■ 

It is clear from (fTol that the probability of rank failure is 
dependent on the number of users L and the finite field size 
p whereas probability of error at the relay (UTi is dependent 
only on the additive noise. 
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VI. Performance Analysis 

In this section, we present the simulations to illustrate the 
performance of the proposed encoding and decoding func- 
tions in terms of (i) the probability of error at the relay, 
which measures error in detecting linear combinations, (ii) 
the probability of error at the destination, which measures 
the probability of incorrect detection of original signals. We 
consider L = 2 users sending out signals to the destination via 
relay. We study the performance of our scheme using different 
residue classes Q v and their corresponding finite fields F p . 
These classes have been listed in Table I giving the residue 
class, corresponding fields and the u and v values to design 
the isomorphism ip in dT)-®. Further, we consider uniformly 
distributed channel gains between all the nodes. For each 
residue class, we make L x 10 4 transmissions from source to 
destination and the decoding of original signals is done after 
every L transmissions. 

Figure 3 shows the comparison of probability of error with 
varying SNR. It can be seen that a higher order finite field (or 
a higher order Q^) gives a higher probability of error at the 
relay for the same SNR. This happens because the source of 
error at the relay is only the additive noise. The impact of this 
additive noise is determined by packing and a higher order 
field will have a denser packing as compared to lower order 
field for same SNR. 

However, at the destination, the probability of error de- 
creases with increasing SNR up to a certain point and then 
it attains a constant value. This is because the overall error 
is contributed not only by the additive noise at the relay but 
also due to the probability of rank failure at the destination. 
The probability of rank failure is independent of SNR ([16} 
and is fixed for any given field size and number of users. The 
probability of error at the destination decreases with increasing 
SNR only up to the point when it becomes comparable to the 
probability of rank failure for a given field size. After this 
point, the error at the relay becomes negligible as compared 
to error due to rank failure and therefore, error probability 
at the destination becomes a constant equal to rank failure 
probability. A higher order partition gives a lower probability 

'Since noise has Gaussian distribution, Pr = Pr(|| z ||> i) = 1 — 



Fig. 3. Probability of error at the relay and at the destination. In all cases, 
L = 2 users are considered. Three different partitions are plotted. 
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e 2<y' A du ) , and the result follows. 



TABLE I 

Finite fields p (where p = 1 mod 4), it (where p = nir* ) and the 

VALUES OF U, V (WHERE UTT + V1T* = 1) 



of error at the destination at high SNR due to lower probability 
of rank failure as compared to lower order partition like 
Also, note that the theoretical union bound estimate given in 
Theorem 1 is reasonably tight. 

VII. Conclusions 

In this letter, we have introduced a concrete scheme to 
implement Compute and Forward relaying protocol using finite 
size signal constellations. We have designed encoding and 
decoding functions using residue class of Gaussian integers 
and used their natural properties of isomorphism with finite 
fields to obtain mapping between signal space and message 
space. We have obtained an analytical union bound estimate of 
probability of error and validated it via simulations. We proved 
that at high SNR, full rank requirement of the coefficient 
matrix plays the key role in determining the end to end 
performance of CF. 
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